Astronomy & Astrophysics manuscript no. RAVE 


©ESO 2012 


September 5, 2012 





Thin disk kinematics from RAVE and the solar motion 

S. Pasetto 1 - 2 , E.K. Grebel 2 , T. Zwitter 3 - 4 , C. Chiosi 5 , G. Bertelli 6 , O. Bienayme 7 , G. Seabroke 1 , J. Bland-Hawthorn 8 , C 
Boeche 2 , B.K. Gibson 9 - 10 , G. Gilmore 11 , U. Munari 6 , J.F. Navarro 12 , Q. Parker 13 - 14 - 15 , W. Reid 13 - 14 , A. Silviero 5 - 16 , M. 

Steinmetz 16 



(N 

o 

(N 

D 
CO 

m 



< 

o 

6 



> 
o 

q 

OV 

o: 

(N 



University College London, Department of Space & Climate Physics, Milliard Space Science Laboratory, Holmbury St. Mary, 
Dorking Surrey RH5 6NT, United Kingdom 

Astronomisches Rechen-Institut, Zentrum fur Astronomie der Universitat Heidelberg, Monchhofstr. 12-14, 69120 Heidelberg, 
Germany 

University of Ljubljana, Faculty of Mathematics and Physics, 1000, Ljubljana, Slovenia 
Center of Excellence SPACE-SI, Askerceva cesta 12, 1000, Ljubljana, Slovenia 

Department of Physics and Astronomy "G. Galilei", Padova University, Vicolo dell'Osservatorio 3, 35122 Padova, Italy 

INAF - Padova Astronomical Observatory, Vicolo dell'Osservatorio 5, 35122 Padova, Italy 

Observatoire astronomique de Strasbourg 11 rue de l'Universit, 67000 Strasbourg, France 

Sydney Institute for Astronomy, University of Sydney, NSW 2006, Australia 

Jeremiah Horrocks Institute, University of Central Lancashire, Preston, PR1 2HE, UK 

Monash Centre for Astrophysics, Monash University, Clayton, 3800, Australia 

Institute of Astronomy, Cambridge University, Madingley Road, Cambridge CB3 0HA 

University of Victoria, Department of Physics and Astronomy, Victoria, BC Canada V8P 5C2 

Department of Physics and Astronomy, Macquarie University, NSW 2109, Australia 

Macquarie research centre in Astronomy, Astrophysics and Astrophotonics, Macquarie University NSW 2109, Australia 
Australian Astronomical Observatory, PO Box 296, Epping, NSW 2121, Australia 
Leibniz-Institut fur Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany 



Accepted for pubblication on A&A 



ABSTRACT 



Aims. We study the Milky Way thin disk with the Radial Velocity Experiment (RAVE) survey. We consider the thin and thick disks 
as different Galactic components and present a technique to statistically disentangle the two populations. Then we focus our attention 
on the thin disk component. 

Methods. We disentangle the thin disk component from a mixture of the thin and thick disks using a data set providing radial velocities, 
proper motions, and photometrically determined distances. 

Results. We present the trend of the velocity dispersions in the thin disk component of the Milky Way (MW) in the radial direction 
above and below the Galactic plane using data from the RAdial Velocity Experiment (RAVE). The selected sample is a limited 
subsample from the entire RAVE catalogue, roughly mapping up to 500 pc above and below the Galactic plane, a few degrees 
in azimuthal direction and covering a radial extension of 2.0 kpc around the solar position. The solar motion relative to the local 
standard of rest is also re-determined with the isolated thin disk component. Major results are the trend of the velocity mean and 
dispersion in the radial and vertical direction. In addition the azimuthal components of the solar motion relative to the local standard 
of rest and the velocity dispersion are discussed. 



Key words. Stellar kinematics and dynamics - 
evolution, thick disk 
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1. Introduction 

The study of disk galaxies has played a fundamental role in 
the attempt to understand the intricate mechanisms that g overn 
galaxy evolution (see, e.g.. Ivan der Kruit & Freemanll201 ll for a 
recent review). The complex challenge of their origin and evo- 
lution has occupied th e interests of t he astronomical community 
for a long time (e.g . . iBertinl l2000fc IBinnev & Merrifield|f l998; 
Comb eset al.l 1 1995b IBinnev & TremaineT 119871: ISaslawlll985l) . 
However, the global picture of the origin and evolution of spi- 
ral galaxies requires thorough comprehension of the structure 
and evolution of our own Galaxy, where many of the relevant 
properties and processes can be studied in detail. In this process, 
by mapping the Galactic disk distribution, obtaining a detailed 
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description of our Galactic disk kinematics is of paramount im- 
portance in order to determine the structure and the evolution of 
the Milky Way. 

The first fundamental steps in this process were done 
by the determination of the characteristic parameters of the 
disk through the study of the star-count equation (e.g., 
Trumrjler & Weaver] 1 19531) as in the pioneer i ng m odels of 
Bahc all & Soneiral (1 19801) or iGilmore & Reidl (Il983l) . These 
studies opened the way for more complete investigations by 
takin g into account also the kinematic compone nt of the data 
(e.g.. iBienavme et alJll987t iMendez et ail 1 19961) . As a conse- 
quence of the constantly growing amount of available data 
and our improving knowledge of Galactic structure, the initial 
ideas of galaxy formation (e.g.. | Eggen et al.l 1962: Sea rle & Zinnl 
119781 : iPartridge & Peebles! [l967) have developed dramatically 
and merged with cosmological research becoming now either 
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a part of a standard textbook picture (e.g., Bin nev & Merrifieldl 
ll998tlCarroll & Ostliell2006l) or a cuiTent hot topic of research. 

Continuing this long line of research the RA dial Velocity 
Exper i ment project (RAVE) was designed (see IZwitter et alJ 
2008; SteinmetzetaJJ [2006): a stellar survey that provides 
spectroscopic parameters for stars of the Galaxy that are 
mainly part of the thin disk, thick disk or halo component. 
The power of this survey has been repeatedly demonstrated 
with a series of papers focused on the principal target of 



matics (e.g. 


Veltz et al. 


2008; 


Casetti-Dinescu et al. 


2011; 


Siebert et al. 


2008; Hahnetal. 2011; Coskunoglu et al 


2011; 


Siebert et al. 


120111 IWilson et al. 


2011; Karatas & Klement 


2012), identifying stellar streams (e.g., Klement et al 


2008; 


Seabroke et al. 2008; Williams et al. 201 1; Kiss et al. 201 1) and 



studying the chemistry of Milky Way (MW) components (e.g., 
Ruchti et al l 1201 Ol 1201 ll iFulbright et all 120101; iBoeche etal] 
20lUICoskunogLuetal.ll2012l) . 



In this paper we proceed with the general efforts devoted to 
understanding the properties of the Galactic disk from a kine- 
matic point of view. Ultimately, a large spectroscopic sample of 
disk stars can reveal the characteristic parameters of the structure 
of our Galactic disk, the evolutionary history of the disk popula- 
tions and the physical processes that acted on the stars in order 
to produce the picture that the RAVE data show us today. 

Our task is to attempt to break the parameter degeneracy that 
results from the imprecise (or incomplete) knowledge of the ob- 
served parameters for the catalogue stars. The RAVE survey pro- 
duces directly several phase space parameters. We will directly 
work with: the Galactic coordinates (I, b), the radial velocities 
v r = (v'hei,?)?, where the unitary vector to a generic star is a 
function of the directions f = f (I, b). Within their errors these 
parameters permit us to constrain the true velocity vector once 
the distance is known. All these parameters come, of course, 
with their uncertainties. In particular, with indirect methods the 
distances ivi = rheir of the stars ca n be derived only within a 
given error (e.g.. lZwitteretalJl2010l) . 

Although these quantities, if complemented by accurate dis- 
tances and proper motions would provide the full phase space 
parameters for each star, i.e., the initial conditions for a complete 
N-body integration or equivalently the boundary conditions for 
the Liouville equation, there are several reasons that prevent us 
from following this direct approach: 

- the error bars of the above mentioned determinations leave 
too large a space of possible solutions, which prevents us 
from achieving an acceptable convergence of the results; 

- the proximity of the stars in our sample prevents us from 
really following the trajectory of each star, not permitting us 
to evaluate the exact behaviour of more distant parts of the 
Galaxy. 

Thus, we are in the apparently paradox situation that success- 
ful methods for external galaxies do not yield the same quality of 
results in our much closer G alaxy (see e.g., the results based on 
the full N-body integrations dFuxll 19971) . on Schwarschild meth- 
ods or t heir more recen t evolution, the Made to Measure methods 
e.g., bv lDehnenl d2009l) ) . The Galaxy description we are going to 
utilize here is only of statistical nature, and it is based on the col- 
lisionless Boltzmann approximation to describe the Galaxy as a 
set of two components (we are considering only thin and thick 
disk stars in our description). This will prevent us from easily de- 
riving any consideration of dynamical nature, and we will limit 
ourselves to a simpler kinematic description. In particular, here 
we present a statistically determined description of the thin disk 



component alone that can be achieved, within reasonable errors, 
in the solar neighbourhood. 

Nonetheless, the problem of extending the kinematic param- 
eters estimated here for the solar neighbourhood to the entire 
disk can be justified based on symmetry considerations and on 
the comparison with external spiral galaxies. If we limit our- 
selves to a kinematic description, i.e., a steady state description 
for a one-comp onent disk (namely the th in disk), then the Jeans 
theorem (e.g.. iBinnev & Trem aine 1987) authorizes us to look 
at the sola r vicinity as a parad igm of the Galactic disk kinemat- 
ics (e.g jMav & Binnevlll98 6 |) as long as we do not attempt too 
detailed a description, (e.g., Ouillen et al. 201 1). We will exten- 
sively discuss the limitation of our approach in later sections of 
our paper. 

The structure of this paper is the following: in Sect. [2] we 
review the hypotheses of the methodology adopted and the data 
selection criteria. In the Sect. [3] we present the results and in 
Section 4 our conclusions. 



2. Methodology adopted 

2.1. Target 

In this paper, our target is to present the trend of the first two 
moments of the velocity distribution (the mean and the central 
moment of order 2) for the thin disk component along the merid- 
ional plane ( Q,R,z) of the MW. Using a methodology recently 
developed in Paset to et al.l (I2012I) (hereafter Paper I) we are able 
to disentangle the thin disk and thick disk components from the 
mixture of the two disk populations near the sun. Then, with 
the help of proper motions and distances photometrically deter- 
mined by Zwitter et al. (2010) we can follow the trend of the 
thin disk velocity ellipsoid described in cylindrical coordinates 
(0,R,<p,z). 

2.2. Data selection 

We review here the data selection, while the methodology is fully 
explained in Paper I. We used an internal release of the RAVE 
data catalogue equipped with photometrically determined dis- 
tances (Zwitter et al. 2010), radial velocities and proper motions. 
This provides us with a full phase space description for about 
210,000 stars consisting of a mixture of thin and thick disk and 
halo stars. 

In order to limit our sample only to the thin and thick disk 
components, we select from the RAVE catalogue all the stars 
with log 10 g > 3.5 (with g denoting the surface gravity of the 
star), the near-infra red colours from the Two Micron All Sky 
Survey (2MASS), ( Skr utskie et all l2006h ./ - K € r0.3 ,l.ll 
mag in the J and K band (see ISteinmetz etal] d2006l) and 
Zwit ter et alJ d2008l) for de tails), the latitude of the sample \b\ > 
10 deg, [Fe/H] > -1 (from lZwitter et a l. (2008) in order to min- 
imize any halo contamination), ||Vh e ill < 150kms~' as absolute 
value of the true velocity vector and S /N > 20 for the signal- 
to-noise ratio of the spectra. Finally, we cut the sample with the 
kinematical condition investigated in Paper I: Av^ AVE (/,-, bi) > 



2.0(5vf (li, bi) (Eqn. (6) in Paper I) where 6v° (/;, b t ) is the differ- 
ence between a maximum value max (v^ and a minimum value 

min(v^) of the Galactic component of the radial velocity in the 
direction of the star s, where re [0, 1] kpc. This condition has 
been investigated in Paper I in order to reduce the indetermina- 
tion of the kinematical parameters for the most distant stars and 
the large scale effects along the selection process. An isocon- 



2 



S. Pasetto et al.: Thin disk with RAVE 




tons" 1 2 4 6 8 10 12 14 16 



Fig. 1. Celestial sphere distribution in an Aitoff projection of the 
selected stellar sample (yellow dots) as outlined in the text. The 
background density map is explained in the text. 



tour plot of this inequality ba sed on the assumed Galactic grav- 
itation al potential field (e.g., Ber telli et al.ll2.Q03t IVallenari et al.l 
2006) is given in the background contour plot in Fig.Q]and can 
be directly compared with Fig. 1 of Paper I. As we can see from 
this contour plot, the highest <5vp (/;, £>,) (yellowish zones) oc- 
cur at low latitude and the smallest in the red zones. For b — 
Svy (/,-, bd is the terminal velocity relation and we see that the 
the largest induced error is of the order of 16 km s _1 at a dis- 
tance of 1 kpc and at low Galactic latitudes not covered by the 
RAVE survey. Here we assume that the true velocity vector is 
defined in a suitable reference system for the velocity space 
(O, U, V, W) aligned with the configuration reference system at 
the solar position. The U direction points to the Galactic coor- 
dinates (I, b) = (180, 0) deg, V runs through (/, b) = (90°, 0) deg 
and W points to the North Galactic Pole (NGP). 

This data selection leaves us with a sample of 38,805 dwarf 
stars that map the solar neighbourhood's stellar distribution 
avoiding the influence from the small halo component (see 
Figure [TJ. 

2.3. Disentangling technique 

In the classical kinematic description, a single collisionless stel- 
lar population satisfies the collisionless Boltzmann equation, i.e., 
a linear and homogeneous partial differential equation in the 
phase space distribution function (DF), / (once the total poten- 
tial (Dtot of the whole Galaxy is assumed). Since this equation is 
linear in /, the whole system can be described as a mixture of 
populations, e.g., / = Yifi- in Paper I (to which we refer the 

i 

reader for further details) a methodology based on the Singular 
Value Decomposition solution for overdetermined systems was 
outlined in order to statistically derive the first moments of a 
DF / m ; x = /thin + /hick for a stellar mixture of thin and thick 
disks. Then, we statistically disentangled the thick disk popula- 
tion alone, assuming a generali zed ellipsoidal DF, /thick, in mod- 
ification of the original work of lCubarsi & Alcobel(l2004|) . 

The selection procedure works only on the radial velocity. 
At every iteration, say j, by randomly choosing a radial velocity 
value within its error bar, say v; (si) for each star s, e P^\ i.e. a 
star that passed our selection criteria in Section |Z21 we are left 
with a selected RAVE sample of stars within their associated 



errors: 

/>« = { S .| V W e [Vr (sj) _ AVr , Vr (si) + AVr j (i) 

The procedure then isolates the thick disk component by ran- 
domly excluding stars in order to maximize a statistical weight 
param eter q described in Paper I (see also ICubarsi & Alcobel 
2004). Once the population that maximizes q within the machine 
error precision (namely q — oo) has been isolated, we are natu- 
rally left with two samples. The sample S ^, with c P^Vj, 

describing the thick disk and the sample P^\S^' or (read 
"not pW") describing the thin disk, is defined as: 

T Ui = pUl\ S lJ] = [ s .\ s . € pin A s . ( s li] }Vj (2) 

thus once the thick disk sequence of stars for the iteration j is 
removed all the other stars are considered to be thin disk stars. 
We explicitly state that: 

- we do not limit the determination of the errors to a prese- 
lected sample of stars. For instance, a generic star s, of P^ 
can belong to the sample for the iteration number j and 
to the sample for the iteration number j + I. The indices 
j and i are not correlated. 

- the thin disk selection is based only on the radial velocities 
and their errors once the sample P has been chosen. This 
makes use of the good quality errors as determined by the 
RAVE survey in Eqn. 1 . 

We emphasize that the first point is fundamental: we are not 
assigning any probability weight to a single star in the sense of 
being a thin disk or thick disk star, even though we know its full 
phase space description. This assigns much more generality to 
our method and prevents us from biasing the results. 

A schematic view of how the method works is sketched in 
Figure |2] where for a generic sample P^\ the direction to each 
star, f , is represented by grey arrows (solid for the thin disk stars 
and dashed for the thick disk). The true, unknown velocity vec- 
tors for the thick disk stars vS ck , drawn from the DF /hick, are 
represented in light red. The sketch shows the thin disk star ve- 
locity vectors, vr? 1 , with dark red arrows that are slightly greater 
closer to the Galactic centre (GC) and smaller in the radial out- 
ward direction for R > R Q , with R G being the solar position. This 
is an expected trend for the thin disk velocity ellipsoid which is 
supposed to have larger ellipsoidal axes i n the inner zones of 
the G alaxy and smaller ones outside (e.g., Binnev & Merrifield 
1998). The radial projection of each velocity vector is shown 
with red arrows: thinner for the thin disk and thicker for the thick 
disk. We used in the selection process only the radial velocity 
(red arrows) vy = (Vhei, ?) ? provided by RAVE, and thanks to 
the different velocity dispersion of the thin and thick disk com- 
ponents (graphically realized with longer and shorter red arrows) 
we managed to disentangle the thick disk component as in Paper 
I. Once the thick disk stars have been selected (dashed grey ar- 
rows) we call all the other stars "thin disk stars" (the solid grey 
arrows). We determine their true velocity vectors (dark red ar- 
rows) by making use of the proper motion values tabulated in the 
RAVE internal data release and calibrated with photometrically 
determined distances from Zwitte r et al. (2010). In addit i on we 
use the matrix transformation of Johnson & Soderblom (119871) 
corrected to J2000 as in iPasetto et al.l ( 20031) in order to derive 
the velocity dispersion tensor for the thin disk alone. 
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Fig. 2. Sketch of the geometry exploited by the method. The 
observer is located at the Sun's position R Q at the centre. The 
Galactic centre (GC) is at the bottom of the figure. The directions 
on the celestial sphere to each star s, are indicated with dashed 
arrows for thick disk stars (the sample S [;] ) and solid arrows for 
the thin disk stars (the sample T^). The red thick arrows repre- 
sent the radial velocities of the thick disk, whose corresponding 
unknown true velocity vector is in light red. The thin disk radial 
velocities are represented by the thin red arrows, and the corre- 
sponding unknown dark red arrows represent the true velocity 
vectors. Note that the dark red arrows are generally longer in the 
GC direction and smaller in the anti-Galactic centre direction as 
expected from the velocity dispersion trend. 



3. Results 

We present now the analysis of the averaged results for a large 
set of statistical realizations of the sample for the single thin 
disk component (typically j = 1, ..,N with N > 1024) as out- 
lined in the previous section. With the concept of a single com- 
ponent we mean each part of the system that is characterized 
by a specific property, i.e. the chemical elements C, the age t, 
the mass-to-light ratio M/L in some specific band, the kinematic 
parameters, the scale length, etc. If we collectively call S the 
sample of the physical properties that we are not interested in, 
its cardinality (i.e. the number of its elements) can be very high, 
or equal to the cardinality of the continuum. We will average the 
DF over all the properties we do not wish to study so that we will 
refer to the thin (or thick disk) DF simply as / t h; n J funnda. 
Since the mapping of the meridional plane kinematics is our 
principal goal, we will limit ourselves to the kinematic descrip- 
tion of the thin disk DF (through only its first and second velocity 
moment). Thus once the sample of stars is selected, we are 
left with the proper motion and their errors, radial velocities and 



their errors, distances and their errors Vs, € (and the sky 
directions considered without errors). 

3.1. Data binning 

Once the thin disk velocity vectors Vhei have been found, we pro- 
ceed by binning the sample in the meridional plane. We are not 
binning in the <p direction in order to exploit the axisymmetric 
assumption we made when the selection procedure was devel- 
oped as explained in Paper I, i.e. we are summing up all the <f> Si 
for each star s, and the Sun at the azimuthal position <p e = 
in a Galactocentric cylindrical reference system. The choice of 
the bin size does not affect the results we are presenting once 
a sufficiently large bin (say, > 100 pc) is adopted since we are 
interested in the general trend of a smooth function of the inde- 
pendent variables R and z (and not in the description of the gran- 
ularity of the DF). Different statistical approaches can be shown 
to give the same p hysical results, e.g., the statistical estimators 
in lVio et alJ (Il994l) . 

3.2. The mean motion 

In the Tables [1] [2] and [3] the velocity components of the mean 
velocity vector are listed along with their errors for a selected 
sample of bins. A common trend is present in the values of the 
tables, although it is also evident that the large error bars prevent 
us from obtaining a stable fitting function. This means that al- 
though the average values for vr, and v z show a smooth trend, 
we will not be in the position to compute their derivatives <9rv, or 
d z Vi which remain a primary target of forthcoming surveys such 
as Gaia. In general the tables show a lack of data in the Galactic 
centre and anticentre directions as a consequence of the selec- 
tion procedure outlined in Section 12.21 and the Galactic latitude 
selection cuts (\b\ > 10 deg). Nevertheless, when possible, we 
retain values as e.g., vr = 18.8 + 45.2 km s _I in Table Q] for 
the distance range z e ]— 0.5, -0.3] kpc (and similar in Tables 
|2] and O, because their mean values are in agreement with the 
general trend of the analysed property, e.g., vr, expressed by the 
neighbouring bins in the Tables. 

3.2.1. Mean radial motion 

The trend of the mean radial motion Vr is presented in Table 
Q] The column describing the vertical range z e ]— 0. 1,0.1] 
kpc contains the closest and our best-determined values in the 
solar neighbourhood. Assuming the location of the Sun to be 
in the range R e ]8.4, 8.6] kpc we recover a value of vr = 
10.9 + 1.0 km s , in gener al agreement with the re cent inde- 
pendent determination (e.g.. lAumer & Binnevl [2009) but with 
a larger error bar due to the statistical approach adopted. This 
value is in good agreement also with our recent determination 
deduced from the thick disk component alone in Paper I. Also, 
the value of vr seems to be robust against possible large scale 
effect correction (see also, Schonrich 2012). Table Q] provide us 
with a consistent value for the LSR radial motion not only within 
the closest radial bin around the solar location, but also in its 
extended neighbourhood. Table Q] shows no net global motion 
within the error bars, i.e., the values show no net global radial 
motion of the stars with respect to the motion of the Sun. This is 
an indication of the absence of a net radial migration trend with 
respect to the solar motion and is in line with results already ob- 
tained much closer inside and fu rther outside of our range of dis- 
tances in the milestone paper of Lewis & Freeman! (TT989). This 
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Table 1. Trend of the mean radial velocity Vr expressed in [km s _1 ] and the distance from the Galactic plane bins are in [kpc]. 
Empty slots indicate parameter ranges with too few stars. 



vr(R,z) 


<-0.6 


]-0.5,-0.3] 


]-0.3,-0.1] 


]-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 




21.7 ±50.0 






63.1 ±0.0 




7.8 


35.3 ± 67.2 


14.2 ± 12.1 


12.8 ±28.0 




8.2 ±21.3 


17.5 ±20.1 


8.0 


3.7 ± 6.6 


16.4 ±7.5 


9.1 ±4.3 




13.2 ±9.3 


14.4 ±9.6 


8.2 


20.4 ± 24.5 


14.9 ±4.4 


9.9 ±2.6 


10.0 ± 11.9 


13.3 ±4.5 


16.4 ± 8.2 


8.4 


19.2 ± 18.7 


11.4 ±3.1 


10.9 ± 2.7 


10.9 ± 1.0 


13.4 ±2.6 


12.5 ±4.9 


8.6 




10.4 ± 3.1 


10.2 ± 1.6 


10.1 ±2.6 


11.5 ±2.4 


5.7 ± 2.6 


8.8 


6.5 ± 11.6 


6.1 ±3.1 


7.4 ±2.2 


17.2 ±45.3 


6.4 ± 4.0 


8.0 ±6.5 


9.0 


-6.2 ±31.7 


6.7 ±3.8 


10.3 ± 10.9 




9.4 ± 26.0 


4.1 ±3.5 


9.2 




-3.8 ± 10.6 


12.9 ±0.0 








9.2 




22.4 ± 0.0 











Table 2. For completemess we present here the trend of the mean circular velocity expressed in [km s ] and the distance from 
the Galactic plane bins are in [kpc]. Empty slots indicate parameter ranges with too few stars. Refer to the text for the limitations of 
our approach along this direction. 





<-0.6 


]-0.5,-0.3] 


]-0.3,-0.1] 


]-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 




6.6 ± 17.1 










7.8 


-43.1 ±78.5 


-12.5 ±8.3 


12.5 ± 26.2 




2.7 ±7.7 


2.2 ±3.6 


8.0 


-47.2 ±61.9 


-10.8 ±3.9 


2.5 ± 1.6 




5.8 ±4.6 


-8.0 ±3.0 


8.2 


-42.6 ±45.9 


-4.5 ± 1.2 


6.0 ± 1.6 


7.8 ±8.8 


6.1 ±2.9 


-8.6 ±2.8 


8.4 


-36.6 ±41.0 




8.3 ± 1.5 


11.3 ± 1.5 


7.2 ±0.9 


-5.7 ±2.0 


8.6 


-42.6 ± 58.0 


-2.7 ±0.5 


9.9 ±2.6 


13.1 ±2.9 


9.2 ±2.2 


-7.8 ±2.9 


8.8 


-44.7 ± 140.4 


-3.0 ±0.6 


9.5 ±4.3 


9.2 ± 29.6 


7.2 ±5.8 


-11.2 ±8.4 


9.0 


-29.6 ± 103.9 


-5.8 ±4.1 


5.0 ±7.7 




5.5 ± 20.9 


-11.7 ±24.9 


9.2 




-2.6 ±4.3 











Table 3. Trend of the mean vertical velocity v z expressed in [km s ] and the distance bins from the Galactic plane are in [kpc]. 
Empty slots indicate parameter ranges with too few stars. 



v-SR-z) 


<-0.6 


]-0.5.-0.3] 


]-0.3.-0.1] 


]-0.1.0.1] 


]0.1.0.3] 


]0.3.0.5] 


7.6 




18.8 ±45.2 










7.8 




15.6 ± 13.1 


15.1 ±32.7 




22.0 ± 56.6 


10.8 ± 13.3 


8.0 


14.9 ± 20.3 


13.5 ±7.3 


10.3 ± 6.5 




6.8 ±4.5 


10.5 ±5.8 


8.2 


18.5 ±22.1 


10.3 ±4.2 


9.0 ± 2.0 


5.4 ±7.9 


8.8 ±2.3 


10.0 ±4.1 


8.4 


13.4 ± 15.5 


10.7 ±4.3 


8.1 ± 1.2 


7.2 ± 1.3 


9.2 ±2.7 


13.4 ±7.2 


8.6 


16.2 ±23.6 


9.9 ±3.7 


9.4 ± 2.0 


8.1 ±3.1 


10.3 ±3.7 


16.2 ±8.8 


8.8 


5.1 ± 16.5 


8.2 ±3.6 


10.8 ±4.3 


10.0 ±33.3 


12.5 ± 11.8 


17.8 ± 18.0 


9.0 


21.3 ±75.9 


9.6 ±9.4 


12.3 ± 19.5 




5.2 ± 20.7 


11.2 ±26.2 


9.2 




25.8 ±47.1 











means that for a mixed set of stellar populations, the assump- 
tion of a zero net expansion is quite reasonable within the error 
bars even if a small outward trend can be evinced from the plot. 
On the plane z — 0, the radial motion of the stars at the solar 
position remains essentially unchanged from -0.5 kpc to +0.5 
kpc. Similar considerations hold along all radial extensions of 
the RAVE data sample within or outside the solar radial position 
for above and below the plane respectively. 

3.2.2. Mean circular motion 

The values for the mean circular motion v^, are presented in 
Table [2] The solar neighbourhood value for this component is 
vy j0 = 11.3 ± 1.5 km s _1 , w hich is slightly higher tha n the 
classical literature values (e.g., Binnev & Merrifield 1998) but 
agrees with recent val ues such as in Pisku nov et alJ d2006l) and 
Scho nrich et alj d2010t) . Nevertheless, this value is slightly more 
difficult to estimate because of the superposition of several ef- 
fects: the Galactic rotation, the peculiar velocity of the Sun rela- 
tive to the local standard of rest at the solar position (LSR), the 
distribution of velocities of the sample of stars, and finally the 



presence of the thick disk stars. Here we have statistically ex- 
cluded the thick disk's kinematic contribution to this value. The 
relation between the other contributors, for a single component 
thin disk model, can be approximated with the help of the Jeans 
equation as 

I Vc (R, Z)\ = [V 2 LSR - |Sf « + *L) + « + °%) 

™\ dR 5z dR j 

with V% R W as m e velocity of the circular orbit at the radius R, 
p as the stellar population density profile and cry the velocity dis- 
persion tensor elements introduced below. This equation shows 
how part of the kinetic energy of any disk population goes into 
rotation and part in the velocity dispersion, which is necessary 
to stabilize th e disk by achieving a suitable Q parameter (after 
lToomrelll964l) . Nevertheless, there are a few important reasons 
why we cannot easily exploit the previous equation in order to 
disentangle the different contributions to the mean circular ve- 
locity: 

1. This result is an expected shortcoming of the methodology 
adopted to disentangle the thin and thick disk components. 
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As already explained in Paper I, although the collisionless 
Boltzmann equation is linear in the phase space distribution 
function, the velocity dispersion moments become linear 
and only if the mean velocities of the partial distribution 
have the same mean (i.e., in our case, we can correct the thin 
disk and thick disk only in the radial and vertical directions). 
In the azimuthal direction clearly this is not the case (as ev- 
ident from Table [2]), preventing us from using the result on 
the dispersion tensor for the thick disk of Paper I in order to 
correctly recover the value vy >G from Equ. (O and Tables [4] 
0and|7l 

2. While the total potential O tot and the scale length, Iir, which 
enter in the density p (R, z) distribution of the populatio n, can 
be rel ated thanks to the Po isson equation (see, e. g., IPasettol 
120051) or the Appendix of IVallenari etall (120061) . unfortu- 
nately the velocity dispersion errors presented below do not 

, dtr 2 „„ , <9<xjL 

permit us to compute the terms and 

3. The age-velocity dispersion relation implies that the spa- 
tial distribution and kinematics of the stars are strongly de- 
pendent on the age of the stars, i.e., the dispersion tensor 
has a strong dependence on the time cry = cr,-,- (R,z, t) as 
well as on the scale height h z = h 7 (t) (e.g.. IWielenlll974l 
119771) . Hence, in equation (f3]l we can express the time de- 
pendence of v<j by writing (R, z, t) = (r, z, crij (R, z, ?)) 
or more briefly = \<p f cry (R, z, f)J to show that if we as- 
sume <s ^jr then the dominant time dependence in \^ is 
given by cr,j = try (R, z, t). In this case, it is convenient to ex- 
trapolate the value of the motion of the solar neighbourhood 

as = lim (cry (^ , z Q , f)) where {fi , z } then is the 

range {R,z} e]8.2,8.4] x ]-0.1,0.1] kpc of Table 2 and 
is the minimal time to consider a y oung population in statis - 
tical collisionless equilibrium (e.g jAumer & BinnevC zOOS*). 
Nevertheless, by applying such a time binning (on the age of 
the stars) we should limit ourselves to the immediate solar 
neighbourhood, namely only the closest bin {^ ,z Q J to deal 
with a suitable amount of stars per bin, which is not the goal 
of the present paper. Furthermore, if the origin of the age- 
velocity dispersion relation stands in the gravitational po- 
tential fluctuation caused by transient spiral structure (e.g., 
iDe Simone et al.ll2004l) the approximation <K may 
fail and the previous argument for the use of Eqn. ® does 
not hold any more. 

Nevertheless, we present these results in Table [2] to be 
complete. We postpone the study of the trend of in the 
meridional plane to fut ure works and different techniques (e.g., 
M oni Bidin et a 1]|20T1 see also the conclusion in Section 4). 

3.2.3. Mean vertical motion 

Finally, the mean vertical motion is presented in Table [3] The 
value around the solar position is 7.2 +1.3 km s in excel- 
lent agreement with the s tandard adopted literature values (e.g., 
Binn ev & Memfieldlll998l) . No apparent trend is evident within 
the errors, taking V/? and Vz into consideration. 

3.3. Velocity dispersion tensor 

The trend of the velocity ellipsoid along the radial and vertical 
direction of the Galaxy is scarcely known due to the solar loca- 
tion within the Galaxy disk z Q = and the resulting dust absorp- 



ohr [kmj x ] 
3d 34 




Fig. 3. An arbitrarily selected element and range from the cry 
tensor, here ctrr as an example of the three-dimensional error 
box that surrounds every point for two different vertical ranges 
(corresponding to the two upper panels of Fig. [4] for the same 
velocity dispersion tensor). Black lines connect the centres of 
the boxes. 

tion that affects the observations towards the G alactic centre. The 
pioneering work of lLewis & Freemanl (f 1989) remains as solitary 
attempt to observationally address this difficult task and we also 
will limit ourselves only to the close solar neighbourhood by ne- 
glecting the use of giant stars in this work. 

Nevertheless, the RAVE survey has the advantage of let- 
ting us work with an extended number of stars, and to permit 
us to sample the trend of the thin disk velocity ellipsoid in a 
radial range of R e [7.6,9.2] kpc from the Galactic centre, 
an extension that is scar cely sampled in other works such as 
Lewis & Freemanl (119891) . 

We consider now the analysis of the terms cr RR and cr zz . 

We start by looking at Fig. [3] where we show the trend of 
o~rr(R,z) with the Galactocentric distance R. The plot shows 
the trends of the data, but due to the overlapping of the error 
bars and the expected symmetry of the data above and below the 
plane caused by the proximity of the Sun to the ideal plane z = 0, 
the plots with the error bars are generally crowded. We present 
an example of these plots only for illustrative purposes and only 
for the element ctrr in order to give to the reader a feeling of 
the behaviour of the error bars. For each data point, the errors 
are given by a triplet of numbers, { AR, Az, Aery (R, z)J, making 
a visualization hard to realize. Each point lies within a three- 
dimensional error cube whose size is mainly a consequence of 
the error propagation formula (an element by element explic it 
derivation can be found in the Appendix of iFuchs et alJ (|2009)) 
and it depends on the number of stars within the selected spatial 
range (which, in turn, is variable for each realization i* 1 ). 

In order to avoid realizing such a complicated three- 
dimensional plot, in the figures we are presenting for a generic 
element, say ctrr(R,z) in Fig. [5] [6] and|7j we are omitting the 
errors in the z direction. 

Moreover, as a general remark, we can observe how the er- 
ror bars are typically smaller for the data sampled below the 
plane thanks to the better sampling of the RAVE data below 
the Galactic plane. Finally, due to the proximity of the sam- 
ple to the solar location, we prefer to present a linear plot in- 
stead of a logarithmic one ( see, e.g.. iLewis & Freemanl [19891 
Ivan der Kruit & Searld 119811) . The resulting trends for all the 
elements of the velocity dispersion tensor for the thin disk are 
shown in Figure [4] 
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Fig. 4. Plot of all the velocity ellipsoid elements cth in the plane (0,R,z). The axes are in [km s _1 ] for the velocity dispersion (y- 
axes) and in [kpc] for the Galactocentric distances (x-axes). The axis labels are omitted in favour of a small indicator of the element 
we are considering in the lower left part of the plot. In order to visualize different z distance ranges, we use a gradation of increasing 
red for the plot lines from the upper row of plot, z e ]0.3,0.5]kpc down to the bottom row for the range z e ]— 0.5, 0.3] kpc: 
the first row of the upper 6 panels refers to z e ]0. 3, 0.5] kpc, the second to z e ]0.1, 0.3] kpc, the central to the Galactic plane 
z e ]— 0.1, 0.1] kpc. The fourth row includes panels for z e ]-0.1, -0.3] kpc and the fifth one to z e ]-0.3, -0.5] kpc. 



Table 4. Trend of the velocity dispersion element <r RR expressed in [km s ]. The radial distance bins (rows in leftmost column) and 
vertical direction bins (column headers) are in [kpc]. Empty slots indicate regions where not enough stars were sampled. 



(r RR (R,z) 


<-0.6 


]-0.5,-0.3] 


]-0.3,-0.1] 


J-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 




36.6 ± 22.9 










7.8 


45.6 ± 18.5 


35.8 ±7.3 


32.7 ± 20.4 




31.2 ± 25.0 


35.3 ± 10.9 


8.0 


43.6 ± 13.3 


34.8 ± 4.0 


32.6 ±4.6 




31.8 ±5.6 


34.0 ± 5.5 


8.2 


41.3 ± 11.1 


36.2 ±2.9 


32.1 ± 1.5 


28.3 ± 10.9 


32.4 ± 1.8 


36.2 ±4.0 


8.4 


41.9 ± 11.4 


35.4 ±2.1 


31.6 ± 1.0 


27.4 ± 1.1 


31.8 ± 1.2 


34.7 ± 3.8 


8.6 


41.7 ± 13.8 


34.4 ± 2.4 


31.0 ± 1.2 


26.3 ± 1.8 


31.9 ± 1.9 


34.7 ± 5.2 


8.8 


40.6 ±31.6 


34.0 ± 4.0 


32.5 ±3.1 


24.8 ±31.6 


32.6 ±7.9 


33.9 ±9.2 


9.0 


40.9 ± 35.3 


33.8 ±7.7 


31.8 ± 14.5 




31.9 ± 37.7 


32.7 ± 22.9 


9.2 




32.3 ± 17.6 











The trend for cr 
original works of 



is presented in 



Table 
198 



Since the 
(see also 



_ van der Kruit & Searld , . . _ 

Antonucc io-Delogull 19911) . it is commonly assumed that, for a 
self-gravitating thin disk at very small z, a single stellar popula- 
tion is vertically isothermal or only slightly non -isothermal (e.g., 
lAmendt & Cuddefordlll99il iBovv et a"Dl2012l) and has a scale 
height nearly independent of the radius. Under this assumption 
and by considering an exponential dependence of the radial light 
profile ( Freem an! 19701) and a constant mass-to-light ratio, the ra- 
dial dependence of the vertical dispersion can be approximated 

as at. (R, 0) = at (R Q , 0) e . A different approach is based 
on the local constancy of Toomre's stability parameter (iToomrd 
[1961 which in epicyclical approximation leads to a radial de- 

pendence of the form cr\ v (R, 0) = cr 2 RR (RQ,Q)R 2 e (e.g., 
lAmendt & Cuddefordi n~99i). Nevertheless, to first order for the 
solar radial distance R - R bo th the relations resul t in a linear 
fit to the plot of Fig.|7](see also lNeese & Yoss|[l988l) . 



A numerical solution for the velocity dispersion profile 

of the midplane at z = of the MW has been found by 

ICuddeford & Amendtl (Il992l) . Their Eqn. (22) is however de- 
pendent on the gradients of cr zz on the meridional plane or on 
the second derivative of the dispersion tensor. Nevertheless, our 
profiles can be used as a boundary condition for that type of an- 
alytical study. 

From the Tables [4] and [9] we can compute directly the 
anisotropy parameter, which we define for the axisymmetric sys- 
tem as j6 (R, z) = -f- ■ This parameter is related to the first ap- 

°m i— 11— i 

proach mentioned above. A plot of /} is presented in Figures[8]|9] 
andHUl 

In Figure QT| we see the trend of the axis ratio of the ve- 
locity ellipsoid (e.g.. lGilmore et aDll989h along the meridional 
plane. In our approach we are dealing with a composite stellar 
population. Thus we do not necessarily expect to match the the- 
oretical predictions that nevertheless remain based on several as- 
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Table 5. Trend of the velocity dispersion element cr^ expressed in [km s ]. The radial distance bins (rows in leftmost column) and 
vertical direction bins (column headers) are in [kpc]. Empty slots indicate regions where not enough stars were sampled. 



cr R4 ,(R,z) 


]-0.5,-0.3] 


]-0.3,-0.1] 


1-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 


-10.6 + 22.9 










7.8 


-9.8 ±7.3 


-10.2 ±20.4 




-10.1 ±24.3 


-9.7 ± 10.8 


8.0 


-10.5 ±4.0 


-10.2 ±4.6 




-10.0 ±5.6 


-9.9 ± 5.5 


8.2 


-9.6 ±2.9 


-10.9 ± 1.6 


-10.2 ± 10.8 


-10.5 ± 1.9 


-9.6 ±4.0 


8.4 


-10.3 ±2.2 


-8.9 ± 1.0 


-9.8 ± 1.2 


-9.3 ± 1.2 


-9.1 ±3.9 


8.6 


-8.8 ±2.4 


-9.4 ± 1.2 


-9.6 ± 1.9 


-8.4 ± 2.0 


-9.4 ±5.3 


8.8 


-8.8 ±4.0 


-9.7 ±3.2 


-9.7 ±31.6 


-10.0 ± 8.0 


-8.4 ±9.2 


9.0 


-8.7 ±7.8 


-9.6 ± 14.6 




-9.8 ±37.8 


-8.3 ±22.9 


9.2 


-9.6 ± 17.7 











Table 6. Trend of the velocity dispersion element cr Rz expressed in [km s ]. The radial distance bins (rows in leftmost column) and 
vertical direction bins (column headers) are in [kpc]. Empty slots indicate regions where not enough stars were sampled. 



o- Rz (R,z) 


]-0.5,-0.3] 


]-0.3,-0.1] 


]-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 


-22.7 ± 22.9 










7.8 


-17.9 ±7.3 


-14.7 ±20.4 




11.7 ±24.3 


17.8 ± 10.8 


8.0 


-17.8 ±4.0 


-9.1 ±4.6 




11.0 ± 5.6 


16.2 ±5.5 


8.2 


-13.0 ±2.9 


-7.2 ± 1.6 


-4.0 ± 10.8 


8.3 ± 1.9 


12.9 ±4.0 


8.4 


-7.8 ±2.2 


-4.0 ± 1.0 


-3.2 ±2.2 


-2.7 ± 1.2 


9.8 ±3.9 


8.6 


-6.1 ±2.4 


-2.0 ± 1.2 


-4.5 ± 2.9 


-4.8 ±2.0 


2.4 ±5.3 


8.8 


2.8 ± 4.0 


5.8 ±3.2 


-4.8 ±31.6 


-9.9 ± 8.0 


-4.7 ± 9.2 


9.0 


8.1 ±7.8 


9.6 ± 14.6 




-12.5 ± 37.8 


-6.4 ± 22.9 


9.2 


15.4 ± 17.7 











Table 7. Trend of the velocity dispersion element cr^ expressed in [km s ]. The radial distance bins (rows in leftmost column) and 
vertical direction bins (column headers) are in [kpc]. Empty slots indicate regions where not enough stars were sampled. 



o-^(R,z) 


]-0.5,-0.3] 


] -0.3,-0.1] 


1-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 


28.3 ± 22.9 










7.8 


27.2 ±7.3 


24.3 ± 20.4 




24.1 ±24.3 


30.3 ± 10.8 


8.0 


26.8 ± 4.0 


26.1 ±4.6 




24.2 ± 5.6 


27.0 ±5.5 


8.2 


24.8 ± 2.9 


23.0 ± 1.6 


20.3 ± 10.8 


21.4 ± 1.9 


26.5 ± 4.0 


8.4 


23.8 ± 2.2 


21.7 ± 1.0 


20.8 ± 1.2 


20.8 ± 1.2 


23.7 ±3.9 


8.6 


24.8 ± 2.4 


21.2 ± 1.2 


19.6 ± 1.9 


21.1 ±2.0 


24.2 ± 5.3 


8.8 


24.4 ± 4.0 


21.8 ±3.2 


18.1 ±31.6 


20.0 ± 8.0 


23.0 ±9.2 


9.0 


24.7 ± 7.8 


20.5 ± 14.6 




17.3 ± 37.8 


22.9 ± 22.9 


9.2 


24.2 ± 17.7 











sumptions. The ratio ^tjt| + 1 is shown in Figure [TT] In agree- 
ment with our initial hypotheses, once / t hi n = /thin (E, L z , I3) is 
assumed for an axisymmetric and stationary thin disk with an- 
gular momentum L z and the energy E, then this trend has to 
be interpreted as observational evidence of the dependence of 
/thin on a third integral of motion ^3 and not on a deviation from 
the kurtosis values expected from the Schwarzschild distribution 
(see Paper I for further details). As can been seen, in the plane 
e [0.4, 0.5] and outside the plane we are close to con- 
stancy within the error bars. Nevertheless, a small trend can be 
observed in Fig. [TT] that can be interpreted as a suggestion of a 
small decrease of this ratio. 

Finally in Table [9] the trend of cr zz = <r zz (R,z) is reported. 
As is evident, the thin disk distribution of the mixture of main 
sequence stars s eems to be remarkably isothermal (see, also 
iBovv et al.l2012f) . This is an apparent effect due to the mixture of 
the different populations acting to hide the velocity dispersion- 
age relation. The study of this relation is beyond the scope of 
this paper but can be easily achieved without a statistical proce- 
dure to disentangle the thin and thick disk component in order to 
consider a more relevant statistical example (or vice versa: age 
and chemistry can also be taken into consideration to improve 
the disentangling technique, see Sect. [4]). 



The coupling between vertical and radial direction is shown 
with the off diagonal term cr Rz = cr Rz (R,z) in Table [6] This rep- 
resents the tilt of the velocity ellipsoid in the radial direction in 
cylindrical coordinates. By symmetry we expect and we confirm 
that on the plane this tilt is approximately zero; cr Rz (R, 0) = 0. 
The tilt angle here has the remarkable meaning to represent a 
projection of the third integral of motion I3 on the surface sec- 
tion of the meridional plane. Using A = j ^arctan ^ J^^i j j, we 

can obtain (from Table [4] and |9]l a smoothly growing trend for 
the tilt as we get to higher |z|, even though the trend is not per- 
fectly zero at our Sun's position z = Zsun- Nevertheless, within 
the errors, the values are in agreement with the independent de- 
termination presented in IVallenari et alJ d2006l) and the work of 
lAmendt & Cuddefordldl991[) (or w ith the independent, compat- 
ible formulation by Bienavme 2009). 

Tables [5] and [7] report the values for cr^ = cr^ (R, z) and 
o"R0 = cr R (R, z) respectively. The first of these terms, <x w = 
(v^, - v^) 1 ^ 2 represents the second central moments of the az- 
imuthal direction distribution. Despite the bias influencing the 
determination, a smooth trend is still visible, which suggests 
an increasing velocity dispersion toward the Galactic centre as 
well as in the direction of increasing |z| with remarkable symme- 
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Table 8. Trend of the velocity dispersion element cr^ z expressed in [km s ]. The radial distance bins (rows in leftmost column) and 
vertical direction bins (column headers) are in [kpc]. Empty slots indicate regions where not enough stars were sampled. 





J-0.5,-0.3] 


]-0.3,-0.1] 


]-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 


5.1 + 22.9 










7.8 


4.8 + 7.3 


3.3 ± 20.4 




5.2 + 24.3 


4.9+10.8 


8.0 


2.9 ± 4.0 


2.1+4.6 




4.5 + 5.6 


7.3 + 5.5 


8.2 


0.6 + 2.9 


4.3 + 1.6 


3.3 + 10.8 


3.4+1.9 


8.5 + 4.0 


8.4 


-1.6 + 2.2 


3.7 ± 1.0 


3.0+1.2 


4.6+1.2 


9.6 + 3.9 


8.6 


-3.0 + 2.4 


5.8 + 1.2 


5.4+1.9 


6.0 + 2.0 


9.6 + 5.3 


8.8 


-3.6 + 4.0 


5.3 + 3.2 


3.1 ± 31.6 


5.6 + 8.0 


10.8 + 9.2 


9.0 


-4.9 + 7.8 


4.8 ± 14.6 




6.5 ± 37.8 


10.5 ± 22.9 


9.2 


-4.7+ 17.7 











Table 9. Trend of the velocity dispersion element cr zz expressed in [km s ]. The radial distance bins (rows in leftmost column) and 
vertical direction bins (column headers) are in [kpc]. Empty slots indicate regions where not enough stars were sampled. 



<r zz (R,z) 


]-0.5,-0.3] 


]-0.3,-0.1] 


]-0.1,0.1] 


]0.1,0.3] 


]0.3,0.5] 


7.6 


26.2 ± 22.9 










7.8 


25.2 + 7.3 


22.1 + 20.4 




22.1 +24.3 


25.5 ± 10.8 


8.0 


23.1+4.0 


20.1+4.6 




20.8 ± 5.6 


24.5 ± 5.5 


8.2 


20.0 ± 2.9 


18.7 + 2.6 


16.7+10.8 


19.4+ 1.9 


20.3 ± 4.0 


8.4 


19.6 + 2.2 


18.0 + 2.0 


16.3 + 2.2 


17.3+1.2 


19.1+3.9 


8.6 


20.0 ± 2.4 


17.5 + 2.2 


14.7 + 2.9 


17.3 ± 2.0 


18.6 + 5.3 


8.8 


18.6 + 4.0 


17.4 + 3.2 


11.9 + 31.6 


15.7 + 8.0 


18.0 + 9.2 


9.0 


16.3 + 7.8 


15.2+ 14.6 




14.9 + 37.8 


17.0 + 22.9 


9.2 


14.6+17.7 
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Fig. 5. Radial trend of the velocity dispersion element cr RR in 
the vertical interval z 6 ]0.3,0.5]kpc (lighter colour) and z e 
]-0.5, -0.3] kpc (darker colour). These vertical range intervals 
were plotted together in order to illustrate the symmetry above 
and below the disk plane (at z = 0). 
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Fig. 6. Radial trend of the velocity dispersion element cr RR in 
the vertical interval z e ]0.1,0.3]kpc (lighter colour) and z e 
]— 0.3, -0.1] kpc (darker colour). These vertical range intervals 
were plotted together in order to show the disk symmetry above 
and below the plane (at z - 0). 



try above and below the plane. The second term, cr R ^ is null only 
in the case of axisymmetry, although in the literature we found 
considerable misuse of this assumption in the theoretical mod- 
els. Nevertheless this assumption does not hold in the case of our 
Galaxy as seen in Table [5] and a greater effort should be made in 
order to relax the axisymmetry assumption in favour of a more 
permissive rotation symmetry about the z axis for the thin disk 
DF, e.g., /thin (<p) = /thin (<P + n )- Actually, no analytical DF relax- 
ing this approximation is known or applied. This term is strictly 
related to the phenomenon of the vertex deviation, i.e. the angu- 
lar tilt of the velocity ellipsoid with the radial azimuthal direction 

by an angle if/ so that tan (2tf/) — , \ . The problem is strictly 

related to the problem of the streaming motions. Its impact on 



the velocity ellipsoid in the solar neighbourhood has been em- 
phasized by several authors fe.g.. lSoubiran et al. 20031 IPehnenl 
Il998lll999l) Streaming motions can be recognized as local mov- 
ing groups. ISeabroke & Gilmorel (120071) show the presence of 
local moving groups com plicates the ag e-velocit y dispersion re- 
lation . As pointed out by IPehnenl (Il998l) and iDe Simone et al.l 
(2004), they affect the vertex deviation. Up to now, streaming 
motions have been studied o nly in the solar vicin i ty to get a local 
map of the velocity space (Bassi no et aHll986l: iHilton & Bashl 
IT9Q iMavodl 1972b I Woollevll 1 9701: IPehnenl 1 9981) . Nevertheless, 
it is not possible to perform a check with the Hipparcos data for 
the solar neighbourhood because the vast majority of the RAVE 
stars are outside of the range the Hipparcos parallaxes as well 
as the studies on the local vertex deviation. Finally the afore- 
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Fig. 7. Radial trend of the velocity dispersion element ctrr in 
the vertical interval z e ] -0.1, 0.1] kpc. 
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Fig. 8. Trend of the anisotropy parameter ft for the vertical 
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]0.1, 0.3] kpc (darker colour). 



mentioned problems with the first moment of the data distribu- 
tion along the azimuthal direction is expected to influence the 
second central moments in this direction. In the left panel of Fig. 
[12] we show a graphical representation of the angular tilt of the 
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]-0.5, -0.3] kpc (darker colour). 
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velocity ellipsoid with the radial azimuthal direction. The RAVE 
data seem to support the idea that this tilt results as an effect of 
the streaming motions, showing a larger azimuthal tilt for stars 
closer to the plane. 

Finally Table|8]shows the values derived for cr^ z - cr^ z (R, z). 
They are considerably smaller than the previous terms or null 
within the error bars. cr^ z = for a single stellar population tells 
us that the azimuthal and vertical velocities are approximately 
decoupled, or equivalently that the short axis of the velocity el- 
lipsoid remains aligned with the cylindrical coordinate system 
(as expected at z = 0). The same results seem to hold for the 
selected sample of dwarf stars in the thin disk. 



4. Conclusions 

The subject of this paper is the derivation the first two moments 
of the velocity distribution for the thin disk component based 
on stellar spectra for stars in the solar neighbourhood observed 
by RAVE. This was realized by disentangling the thin and thick 
disk; the thick disk has been the subject of our first paper (Pasetto 
et al. 2012, submitted) and in the present paper we exclusively 
focus on the thin disk. Our main goal is the investigation of the 
solar motion relative to the LSR and the thin disk velocity dis- 
persion trend. We derived two of the solar motion components 
as 

[vu,e,v w ,o} = {10.9 ± 1.0, 7.2 ± 1.3} km s" 1 



10 



S. Pasetto et al.: Thin disk with RAVE 



c 

O 
€> 



GC 
GC 
GC 
GC 
GC 



8300 8400 8500 
R[pc] 



600 
400 
200 



& 
N 

-200 



-400 



8 ^00 8400 8600 
Rftjc] 



Fig. 12. Velocity ellipsoids in the solar cylinder as deduced from 
the velocity dispersions of Fig. [4] from two different points of 
view. On the left, we see the radial-azimuthal tilt of the veloc- 
ity ellipsoids with respect to the direction of the Galactic centre 
(GC). The colours of the velocity ellipsoids refer to the height 
above or below the Galactic plane as indicated in the right panel 
(as well as in the vertical sequence of panels in Figj4j». The right 
panel shows the vertical tilt of the velocity ellipsoids above and 
below the plane in the solar cylinder. 
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Fig. 13. Simple representation of the results on the velocity dis- 
persion trend. The velocity dispersion tensor in each spatial point 
(R, z) can be represented with an ellipsoidal surface defined by 
crJ^ViVj = 1 (by summing over repeated indices). Here we plot- 
ted the velocity ellipsoids with their axes tabulated in the pre- 
vious Tables gl □ and their tilt listed in Tables |5] and 
The colour gradient ranges from light colours above the plane to 
darker colours below the plane. 



All the results on the thin disk velocity ellipsoid can be summa- 
rized in a single figure, Fig. Qj] Here we see the radial and ver- 
tical extensions that we have been able to map with the RAVE 
data as well as the ellipsoid axis ratios and their tilts. All the 
values come from the Tables |4] to [9] 

There are several points that need a more extensive investi- 
gation. We just mention a few here: 



- A methodology to obtain constraints on the disk potential 
and the radial velocity sample. A classical approach to this 
problem is based on a pre-assumed form for an analytical 
distribution function based on three isolating integrals of mo- 
tion / t hin = /thin ihih, h) from which by a projection proce- 
dure on the phase space, one can deduce analytically observ- 
able quantities directly linked to the potential. Unfortunately, 
these approaches, as the one adopted in Paper I, are subject 
to assumptions (axisymmetry, stationarity) that fail already 
in the solar neighbourhood: the most obvious cas e is perhaps 
the v ertex d evi ation but see also re cent work by iLuck et all 
(e.g.. 1201 1 ) and Bal ser et al.l d201 lb . Hence it would be use- 
ful to have an analytical model able to relax these assump- 
tions in order to fully exploit the possibilities of the RAVE 
survey. 

- The relation between stellar evolution, chemistry and the 
kinematics (the velocity dispersion) has been neglected in 
this paper. The reason is that in order to obtain a reason- 
able amount of stars and a reasonable precision in the dis- 
persion tensor determination by mapping of the meridional 
plane, we should work with a limited subsample of stars 
closer to us. Nevertheless, it is in principle possible to further 
narrow down the thin disk component with sel ection crite- 
ria bas ed on chemical arguments. For example. iNeves et"ail 
(2009) determined the trend of abundance ratios as a func- 
tion of [Fe/H] from 451 high-resolution spectra of F, G, 
and K main-sequence stars in the solar neighbou rhood con- 
firming a dichot omy in \a/Fe] found already bvlReddv et al.l 
(e.g., 120061) or iFuhrmannl d2004l) . iNavarro et all (1201 ll) in- 
vestigated the possibility of disentangling thin and thick disk 
stars using a combination of [Fe/H], [a/Fe] and the heavy 
element Eu for a sample of 306 stars. Models mixing kine- 
matics and chemistry are currently in development (e.g., 
ISchonrich & Binnevll2009l) . Their results require further in- 
vestigations with respect to the role of the radial migra- 



tion of stars (e.g ., iBrunetti et al. 201 ll iBensbv et all 



iBirdet all 1201 2^ or gas (e.g.. ISpitoni & MatteuccT 
and the presence of metallicity gradients. Nevertheless, the 
RAVE catalogue is surely suitable for this kind of investiga- 
tion since more and more data are being obtained and an ex- 
tension of this disentangling technique to include giant stars 
and chemical constraints is worthy of investigation. 
- Although the average values for v«, and v z present a 
smooth trend, we are not yet able to compute their derivatives 
<9rV; or d z Vj, which remain targets of pri mary importance o f 
forthcoming surveys such as Gaia (e.g., iRobin et"aT1l2012l) . 
These elements are of fundamental importance in order to 
investigate the shape of the thin disk distribution function in 
order to clarify the importance of the third integral of mo- 
tion over the deviation from the expected kurtosis value for 
a Schwarschild DF. 

For the time being, however, the direct comparison of this obser- 
vational work can be done only with other observational works. 
Finally a few cautionary remarks on the limit of our results. As 
already commented in the introduction, the extrapolation of the 
solar neighbourhood kinematic description to the whole Galactic 
disk relies on the validity of the Jeans theorem in the integral 
description of the thin disk distribution function. This is sub- 
ject to the validity of the stationary condition that holds in a 
broad sense for the group of stars sampled by the RAVE sur- 
vey. Nevertheless, we recall that several processes may cause 
small violations of this assumption in the solar neighbourhood, 
for instance: 
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- Spiral structures and the bar redistribute angular momentum 
acros s the disk and heat it up (e.g., ICarlberg & Sellwoodl 
1985) by causing small deviations from the steady-state 
situation and potentially are indeed a main source infiu- 
encing the second order moment of the / t hi n here studied. 
iMinchev & Famaevl d2010h show an example of mixing by 
combined effects of bar and spirals even if the effect that 
they registered seems to be short-lived. 

- The Galaxy accretes angular momentum from the remnants 
of external dwarf galaxies that are merging with it and the 
accretion is expe cted to lead to significant reorientation of 
the sp in axis (e.g. jBinnev & Mavll986tlSellwood & Binnevl 
2002) by potentially influencing the trend of the first moment 

Of /thi„. 

Other processes that are not naturally considered in the station- 
ary description include, e.g., the influence of the Outer Lindblad 
Resonance of the Galactic bar in the solar neighborhood (e.g., 
lDehnen|[l9 99) or the consta nt stripping of small st ellar streams 
from satellite galaxies (e.g., iHelm i & White] IT999I) . These pro- 
cesses are expected to influence the data in the RAVE catalogue 
when also the giant stars are taken into consideration. 
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